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Abstract 

Many manifestly non-Hermitian Hamiltonians (typically, VT— symmetric complex 
anharmonic oscillators) possess a strictly real, "physical" bound-state spectrum. 
This means that they are (quasi-)Hermitian with respect to a suitable non-standard 
metric 0^7. The domain T> of the existence of this metric is studied here for a 
nontrivial though still non-numerical four-parametric "benchmark" matrix model. 
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1 Non- Hermit ian observables with real spectra 



Quite a few realistic quantum models are characterized by a mere "fragile" stability of 
their bound states. For example, the reality of the energies of certain nuclear-physics 
models may be lost after an "inessential" change of its coupling strengths [1]. Such 
a phenomenon is rendered possible when the Hamiltonian merely remains Hermitian 
with respect to a nonstandard scalar product, i.e., with respect to a "nontrivial 
metric" 6 ^ I in an underlying "physical" Hilbert space. This means that all the 
operators of observables H must obey an unusual rule 

= QHQ~ 1 ±H, e = @t >0 . (i) 

Such a property of (typically: Hamiltonian) H guarantees the reality of the spectrum 
and may be called quasi-Hermiticity. 

Practical objections may occur against the use of a nontrivial metric / /. 
One of their sources lies in the fact that many phenomenological models in Quan- 
tum Mechanics are based on a differential- operator realization of their Hamiltonians. 
Thus, whenever one works with a one-dimensional differential-operator Hamiltonian 
(in units Ti = 2m = 1), 

d 2 

H = — — + V(x) =H\ x E (-00,00) (2) 
dx z 

one usually prefers the most economical Q = I scenario since even its simplest alter- 
natives require a number of additional mathematical considerations [2]. Moreover, up 
to the very recent past it has been intuitively expected that the transition to any non- 
Hermitian generalization of the class of differential-operator Hamiltonians (2) would 
be accompanied by a complexification of their spectrum and by a decaying-state 
re-interpretation of the corresponding wave functions. The latter misunderstanding 
even survived the publication of a few studies [3, 4] which paid attention to specific 
examples 

d 2 

H = -— + U(x)+iW(x)^H\ x £ (—00, 00) . (3) 

where the two independent real potentials U(x) = U(—x) and W(x) = — W{— x) 
happened to generate the real bound-state spectra [5] and where their spatial sym- 
metry/antisymmetry has been re-interpreted as a combined parity-plus-time-reversal 
(called '"PT") symmetry of the Hamiltonian [6]. 

The scepticism (well sampled, e.g., on the Streater's webpage [7]) did not fully die 
out even after publication of several analytic, semi-classical and numerical studies 
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of some new VT— symmetric models (3) in 1998 (cf., e.g., [8, 9, 10]). Their au- 
thors demonstrated that the entire spectrum seems to remain real in comparatively 
large domains V of parameters. At present, fortunately, we witness a reconcili- 
ation and final acceptance of the idea that the complex differential Hamiltonians 
may possess the real bound-state spectra, indeed. The two sets of the fresh 2006 
"state-of-the-art" reports may be found in the August dedicated issue of J. Phys. A: 
Math. Gen. (vol. 39, number 32, pp. 9963 - 10261) and in the September dedicated 
issue of Czech. J. Phys. (vol. 56, number 9, pp. 885 - 1064). 

One of the paradoxes accompanying such a development of the subject is that in 
contrast to the popularity of the various anharmonic-type models (based, first of all, 
on their high relevance in field theory [11]), much less attention has been paid to the 
finite-dimensional matrix versions of the non-Hermitian quantum Hamiltonians [12]. 
One of the reasons is that in the matrix models (encountered, quite naturally, in 
variational calculations [1]) one usually deals with "too many" independent matrix 
elements. The selection and/or preference of some of them might look "too arbitrary" 
in the context of physics and/or "too ambiguous" in the language of mathematics. 

In what follows we intend to fill the gap and to study a model which tries to 
circumvent both these "traps" by containing just a few "relevant" free parameters 
and by being still surprisingly rich in its mathematical structure and consequences. 
Moreover, its "derivation" from the differential operator (3) (cf. section 2) gives it 
a certain generic character while leaving it still purely non-numerical and exactly 
solvable. Last but not least, the not quite expected closed-form feasibility of its 
mathematical analysis (cf. sections 3 and 4) is quite well matched by some of its 
appealing phenomenological features, a few remarks on which are added here in our 
final section 5. 

2 Matrix toy models 
2.1 Variational origin 

In order to interconnect the differential and finite-dimensional A r — state Hamiltonians 
let us start from their most elementary differential harmonic-oscillator example H 
with the eigenstates | n), n — 0, 1, 2, . . . or, in the coordinate representation, 

(x\ 2m) = M {m , +) e x2 / 2 H 2m (x 2 ), (x\ 2m + 1) = xU {m ,-) e x " /2 H 2m+1 (x 2 ) (4) 
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where the symbols Tt n denote the well known Hermite polynomials [13]. The sub- 
scripts ± in the normalization factors J\f( m ,±) are added to emphasize that our basis 
states (4) are, simultaneously, eigenstates of the operator of parity V with eigenval- 
ues (— l) n . The action of the complex conjugation T (mimicking the time reversal 
[8]) preserves these basis states once the normalization factors J\f( m ,±) are chosen real. 

The VT— symmetry HVT = VTH of a given Hamiltonian (say, (3)) with real 
spectrum enables us to normalize all the eigenstates | ip n ) of H in such a way that 

VT (x\ lj) n ) = +(x\ V> n > . (5) 

In effect, the fixed parity of our harmonic-oscillator basis (4) is generalized to the 
VT — symmetry of eigenstates. Once we accept such a convention (speaking about 
"unbroken VT — symmetry of wave functions") we may decompose 

|V) = |^+>-i|^-> (6) 
with real expansion coefficients <p m and Xm in the variational ansatz of the form 

N+ N- 

| V+> = E I 2m ) 0m, I V>-> = E I 2m + !> Xm, N + + N_ = N^OO. 

m=0 m=0 

The partitioning of our basis with | 2n) = | n+) and | 2m + 1) = | m_) and its vari- 
ational truncation with N ^> 1 transform Hamiltonian (3) into a finite-dimensional 
partitioned complex matrix 

S \B 
iC L 

where the untilded letters denote the submatrices with real matrix elements, 

d 2 



H 



Smn — ( 2m 



2m + 1 



dx 2 
d 2 



+ U(x) 
+ U(x) 



2n) = (s T ) , (7) 
2n+l) = (L T ) , (8) 

/ V / mn 



dx 2 

B rnn = (2m| W{x) \2n+l) ^ (B T ) = C mn = (2m + 1| W{x) \ 2n) . (9) 
The superscript T denotes transposition. 
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2.2 Real matrix Schrodinger equations 

After we insert (6) in Schrodinger equation H\ip) = E\ip) we reveal that the re- 
sulting partitioned matrix form of Schrodinger equation for bound states is real and 
non-Hermitian, 




No mathematical contradiction appears in the latter picture since the metric O nat- 
urally becomes singular on the boundary dT> of the domain T>. 

In general, the solutions of eq. (10) must be constructed numerically. The well 
known exception is represented by the two-state models [14]. The two-dimensional 
version of our present simplified eq. (10) has also thoroughly been discussed in our 
recent letter [15]. In the corresponding two by two matrix Schrodinger equation 

(:;)(:)-(:) 

parameters s, b and / are real and, by assumption, the even and odd unperturbed 
energies are nondegenerate, s ^ I. Via a suitable scaling we may achieve that 
I — s — 2. A shift of the energy scale E — > E + const leads to the completely 
symmetric arrangement of our entirely general H with s = — 1 and 1 = 1. We recall 
the secular equation det(H — E) = and deduce the energy levels, 

E = E± = ±Vl - b 2 . (12) 

Thus, for our single-parametric family H = H(b) of the N = 2 Hamiltonians the 
domain T> = D(N) of the (single) free parameter b where the energies are real 
coincides with the (closed) interval of b G T>(2) = [— 1, 1]. At both the ends of this 
interval our Hamiltonian ceases to be diagonalizable. For this reason the domain of 
the quasi-Hermiticity of H is often being re-defined as a mere open set or interior 
p(°)(2) = (—1, 1). Under both these conventions, one finds complex energies in the 
vicinity of every element of the boundary dV = dV^ [16]. 

2.3 Anharmonic-oscillator-like four-by-four matrix model 

In the harmonic-oscillator model itself the evaluation of the matrix elements remains 
trivial and one arrives at the simplest illustrative example H containing just a 
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decoupled pair of diagonal submatrices, 



S£i = (2m | H | 2n) = 5 mn ■ {An + 1), (13) 

L^ n = (2m + 1 1 H | 2n + 1) = 5 mn • (4ra + 3). (14) 

For all the Hermitian generalizations of Hq with unbroken parity (V H = HV) equa- 
tion (10) would stay decoupled (B = B T = 0). This means that the parity-preserving 
and Hermitian anharmonicities may be considered "trivial" in leaving the matrices 
S and L decoupled and diagonalizable by the separate unitary transformations in 
the respective even-parity and odd-parity subspaces. 

We intend to employ just the diagonalized and purely harmonic submatrices 
(13) and (14), studying merely the role of the off-diagonal anharmonic-oscillator- 
like coupling matrices B in what follows. Thus, we shall start from the general 
VT— symmetric model (10) with 



H 










B12 ■ ■ ■ 





5 


B21 




-B n 


—B 2 i . . . 


3 


... 


—B\2 
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\ 



(15) 
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We shall restrict our attention to the "first nontrivial" four-level system with the 
truncated dimensions iV + = iV_ = 2. For the sake of convenience we shall also 
symmetrize the unperturbed spectrum via a shift of the origin on the energy scale, 
(1, 3, 5, 7) — > (—3, —1, 1, 3) and arrive at the Schrodinger-equation 



V 



-3 





c 


b \ 


/ 00 \ 





1 


a 


d 


01 


— c 


—a 


-1 





Xo 


-b 


-d 





3/ 


V xi J 



= E 



Xo 

V xi ) 



(16) 



In comparison with the current two-state analyses, a combined effect mediated by 
the simultaneous growth of all the four real parameters a, b, c and d will be more 
complicated of course. At the same time, the levels coupled by an off-diagonal matrix 
element will still follow the pattern revealed at N — 2. This means that, say, the 
growth of c will cause a mutual attraction of the energy levels —3 and —1, etc. 
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Obviously, the separate effects of attraction will compete. One even might encounter 
the usual crossing of levels, not accompanied by any instability and/or subsequent 
complexification of the pairs of the levels involved. This is the reason why the "first 
nontrivial" N = 4 model deserves a deeper analysis. 



det 



0. 



(17) 



3 Constructive analysis of the four-by-four model 

A priori we may say that the influence of the variation of all the quadruplet of 
coupling constants in (16) is tractable non-numerically since the spectrum of energies 
coincides with the set of roots of the secular determinant 

( -Z-E c b \ 

1 - E a d 

-c -a -l-E 

\ -b -d 3-E J 

The exact energies remain obtainable using closed formulae since the corresponding 
secular polynomial is of the mere fourth order, 

E 4 - (10 - a 2 - b 2 - c 2 - d 2 ) E 2 -A (c 2 - d 2 ) E + C(a, 6, c, s) = (18) 

where we abbreviated 

C(a, 6, c, d) = 9 - 9 a 2 - b 2 + 3 c 2 + 3 d 2 + a 2 b 2 + c 2 d 2 - 2 abed . 

Still, the use of the closed formulae does not facilitate our insight in the structure of 
the spectrum too much as it proves prohibitively uncomfortable. Our experience is 
that virtually any alternative analytic approach to eq. (18) proves preferable. 



3.1 Quadruple mergers of the energy levels 

We intend to describe the mechanism of a complexification of the energies in the 
manner which would separate the essential and inessential influence of the variations 
of the parameters. Thus, in a formal language we shall search for the values of the 
matrix elements a bis d at which an abrupt, qualitative change of the spectrum could 
occur. 

In this sense, the most interesting situation occurs at the "points of maximal 
nonhermiticity"'(PMN) at which all the four energy levels coincide, E — E ± — E 2 — 
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E 3 = z = z( PMN \ In the light of an "up-down" symmetry of our unperturbed 
spectrum (—3, —1, 1,3) we may fix 2:( PMA = 0. A change of this value could only 
be caused by a (presumably, perturbative) modification of our model. 

Under the assumption z = our secular equation should read (E — z) 4 = E 4 = 
so that the quadratic term in eq. (18) must vanish, 

a 2 + c 2 + b 2 + d 2 = 10. (19) 

This means that all the four PMN parameters must lie on a four-dimensional sphere 
with radius y/lO. Similarly, from the condition of the vanishing of the linear term 
we deduce that c 2 = d 2 . Finally, the condition C(a, b, c, d) — reads 

9 - b 2 - 9 a 2 + 3 d 2 + 3 c 2 + c 2 d 2 - 2 cdba + a 2 b 2 = 

and degenerates to the factorized relation 

C{a, b, c, d) = (d 2 -ab + zf -(b- 3a) 2 = (d 2 - a) (d 2 - (5) = (20) 

where a — (b + 3) (a — 1) and (3 = (b — 3) (a + 1). This means that at any fixed value 
of d 2 > we get all its solutions (a, 6) as points in the a — b plane which lie on the 
four branches of the two hyperbolas d 2 = a(a, b) and d 2 = /3(a, b) as displayed for 
illustration in Figure 1, with their two centers marked by the bigger circles and with 
the two intersections marked by the small circles (units and axes are dropped here 
as irrelevant). 

Once we return to the former constraint (19) we may conclude that the points on 
the hyperbolas are spurious unless they lie also on the centered circle with the radius 
VlO — 2 d 2 . Hence, under our spectrum-symmetry assumption z = there exist four 
PMN matrix-element solutions which induce a "maximal" , quadruple merger of the 
real energy levels, in a finite interval of values of the free parameter d 2 of course. 
In our illustrative Figure 1 (where we choose d 2 = 8/5) we see that and how the 
resulting points of the boundary dV in the a — b plane (denoted by symbols C2a, 
C2b, C5a and C56) emerge as intersections of the central circle with the respective 
hyperbolic segments C2 — C3 and C5 — C6. 

3.2 Simplified four-by-four model with c 2 = d 2 

In terms of the abbreviations 

A = 5 - d 2 - - (a 2 + b 2 ) , B = (d 2 - ab + 3) 2 - (b - 3 a) 2 
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the symmetry assumption c 2 = d 2 makes our original secular eq. (18) simpler, 

E 4 — 2AE 2 + B = 0, 
and much more easily solvable by the compact formula, 



E±,± = ±\/A±VA 2 -B. (21) 

This means that the necessary and sufficient condition of the reality of the energies 
is given by the pair of requirements 

A>0 (22) 

and 

A 2 > B > . (23) 

Conditions (22) and (23) represent an exceptionally transparent implicit definition 
of the quasi- Hermiticity domain V and/or of its boundary set dV of all the complex- 
ification points. 

Complementing the discussion presented in paragraph 3.1 above we might notice 
that B = C(a,b,d,d) in our older notation. This means that the two hyperbolas 
of Figure 1 represent precisely the boundary curves of the domain D of validity 
of condition B > of eq. (23). One can easily verify that its subdomain where 
d 2 < a(a,b) and d 2 < (3(a,b) consists of two disjoint subsubdomains D(+, A/ B) 
with the respective boundary curves Al — A2 — A3 and Bl — B2 — B3. Similarly, the 
second, single and simply connected subdomain D(— , C) of D where d 2 > a(a, b) and 
d 2 > j3(a, b) is specified by its two pieces of boundary C1 — C2 — C3 and C4 — C5 — C6 
in Figure 1. Obviously, just the latter subdomain has a non-vanishing overlap with 
the interior of the circumscribed circle (19). 

We can summarize that the bound-state energies of the model can only remain 
real inside the latter overlap. In order to arrive at a corresponding sufficient condi- 
tion, one has to recall the last constraint A 2 > B of eq. (23). In its entirely explicit 
form it reads 

(8 + a 2 -6 2 ) 2 >4d 2 [l6-(a + 6) 2 " . 

Its exhaustive discussion and geometric interpretation gets facilitated and becomes 
more or less elementary in its alternative representation 

(2 + a5) 2 >d 2 {A- a 2 ) . 

in the new, rotated coordinates a = a + b and 5 = a — b. 
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4 Special case: VT— symmetric band matrices 



4.1 Perturbative considerations 



H(a, b, c, d) 



(24) 



A re-numbering of the basis (i.e., an interchange of its second and third element) 
makes the matrix in eq. (17) equivalent (i.e., isospectral) to another Hamiltonian, 

/-3 c b \ 
-c -1 -a 
a 1 d 
\ -b -d 3 J 

Once the coupling of the most distant levels vanishes, 6 = 0, and once we re-install 
the symmetry c = d, we arrive at a perceivably simpler two-parametric Hamiltonian 

/-3c \ 
— c —1 —a 
a 1 c 
\ -c 3 J 

It is particularly suitable for perturbative analysis. For example, its one-parametric 
special case 



H(a,c) 



(25) 



H(a) = 



-3 2a 
-2a -1 
-2a 



\ 

2a 

1 2a 



V -2a 3 y 



(26) 



possesses the easily evaluated energies 



E ±1 = ± [-6a 2 + 5- 2 (5a 4 - 12a 2 + 4) 1/2 

£ ±3 = ± [-6 a 2 + 5 + 2 (5 a 4 - 12 a 2 + 4) 1/2 
In the regime of a small a 2 the quickly decreasing curve 



1 1/2 
1/2 



|£ ±3 | = 3_2a 2 -a 4 - -a 6 + a 

6 v 

gets closer and closer to the slowly increasing curve 

\E ±1 \ = l + a 4 + ^a 6 + 0(a 8 ) . 
The energy curves finally intersect, pairwise, at a certain critical strength, 



a 



(CS) 



E. 



(CS) 



E. 



(CS) 
±3 



±1 



±1.612451550. 



Beyond this boundary, i. e., at a 2 > 2/5, all the four energies become complex. 
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4.2 Facilitation of the construction of the metric G 

There exists a clear contrast between the robust reality of the energies resulting 
from a Hermitian Hamiltonian H = W and the globally fragile character of the 
reality of the spectrum in the models which are non-Hermitian and, in particular, 
VT— symmetric. We emphasized in section 1 that this contrast finds a formal rep- 
resentation in the transition to a nontrivial physical metric 0^7. 

On the formal level the operator G may be different for different Hamiltonians 
so that both the Hamiltonian H and the metric G may depend on certain variable 
parameters. One expects, in particular, that the spectrum of H ceases to be real 
out of the domain V of these parameters. Of course, a necessary deeper study of all 
these possibilities is much easier at finite dimensions N when the linear equation (1) 
determines all the eligible metrics G. 

The straightforward linear-algebraic construction of G remains ambiguous. For 
our present, drastically simplified N = 4 input Hamiltonians H the complete solution 
and discussion of the problem remains feasible. For illustration let us consider the 
one-parametric model (26) and solve the related problem 

H*(a)G = GH(a) (27) 

by brute force. This gives the following nontrivial four-parametric real symmetric 
matrix solution 

I ©n ©12 r © 14 \ 

©12 V ©23 S 
r ©32 ©33 ©34 
\ ©14 S ©34 q J 

of the sixteen quasi- Hermiticity conditions (27). In the solution which is routine we 
may employ the notation 

©ii = l(-9p + 3q+l0r + s) + ^(2r-s), 
o cr 



Q(p,q,r,s) 



(28) 



©14 = -t±^h, 33 = 1 (_ 3p - 3 q + 4 r + s) + ± 

6 a 2 



a 1 
©12 = -{3p-3q-4r-s) + -(-2r + s), 
o a 

cx s 
©23 = 7r(-3p + 3g + 2r-s)--, 

a 

01 S 

O u = -(3p-3q-4r-s)--. 
o a 

which specifies the unindexed matrix elements as independent parameters. 
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4.3 Construction of the surface dV near the nonperturbative 
PMN regime 

In the light of our previous results, Hamiltonian H(a, c) of eq. (25) possesses the 
quadruply degenerate energy E = E^ PMN ^ = at the four PMN points with co- 
ordinates a = a (PMiV) = ±2 and c = C (PM7V) = ±y/S. In the vicinity of one of 
them (let us pick up, say, the lower left one) we may set a = a( PMAr ) (—1 + a') and 
c = C ( PMA (—1 -|- c') with some small measures of deviation a' and d . 

In the zeroth order of perturbative analysis this ansatz just reproduces the PMN 
solution a' — c' — 0. On the first-order level of precision the result a' = c' remains 
indeterminate. We have to switch to an improved ansatz containing a new, auxiliary 
small parameter t, 



l+t + at 2 + 0(t 3 )} , c = c( PM ^ [-l + t + 1 t 2 + 0(t 3 ) 



Its insertion in the polynomial secular equation det[H(a, c) — E] =0 (which is of the 
second order in s = E 2 ) leads just to a re-arranged version of the solutions derived 
in paragraph 3.2. In particular, on the second order level of precision we obtain the 
following simplified version of eq. (22), 



10 1 + (-5 + 4 a + 6 7) t 2 + O (t 3 ) > 



which only requires that our small parameter must be non-negative, t > 0. The 
second half (23) of the implicit definition of the quasi-Hermiticity subdomain V in 
the a — c plane is more informative and gives the final, comprehensive estimate 

1 + 1 + 0(t) > a > 7 - \ + 0(t) . 

This formula characterizes the "allowed" parameters a and c which remain compat- 
ible with the reality of the energies. Its form is suitable for the parametric graphical 
plotting of the boundary dT>. The result is sampled in Figure 2 showing that in the 
vicinity of the PMN matrix elements the domain T> has the shape of an extremely 
narrow spike. Its vertex (a( PMN \ c^ PMN ^ represents the simultaneous maximum of 
the size of these elements, saturating the circumscribed-sphere inequality (22) at the 
same time. 
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5 Towards more- dimensional models 



A broad class of modifications of the standard harmonic oscillators may be charac- 
terized by a certain user-friendliness of their perturbative study. A priori, this expe- 
rience may be extended to the quasi-Hermitian models where their N— state matrix 
Hamiltonian is just a small perturbation of the ordinary harmonic oscillator. Beyond 
this perturbative regime, unfortunately, the effects of the non-Hermitian components 
become less predictable. Firstly, in contrast to the usual textbook quantum theory 
where 6 = J, our present use of 6 ^ I (i.e., of a manifest non-Hermiticity of H) 
may mean that the domain V (where H represents an observable) is finite and that 
many of the textbook perturbation-theory theorems and algorithms may cease to be 
applicable [16]. 

In particular, our present study of a specific four-state toy model revealed that 
certain deeply non-perturbative mathematical as well as physical phenomena may 
occur along the boundary &D. Thus, we may expect that perturbation theory can 
offer a reliable qualitative description of the most relevant consequences of the vari- 
ation of the matrix elements only in the regime far from the boundary dT>. In its 
vicinity, on the contrary, perturbative considerations must be used with much more 
care and in an accordingly modified form. 

Several purely theoretical questions emerge near dT> also in the areas of non- 
quantum physics exemplified, say, by magnetohydrodynamics [17], cosmology [18], 
crystal optics [19] or statistical physics [20]. In parallel, the points of &D play an 
important role in the purely mathematical framework of perturbation theory [21] or 
supersymmetric considerations [22]. For all these reasons our present constructive 
study of the boundaries &D may prove relevant in many different applications, after 
an appropriate generalization of our schematic model if necessary. 

In this context, our study of the first nontrivial N = 4 model offered several 
useful hints. We saw that our understanding and reconstruction of the shape of the 
boundary &D will play a key role in the appropriate necessary modifications and 
applications of perturbation techniques. In such a context, it is of course unpleasant 
that the number of the relevant matrix elements (i.e., of the freely variable parameters 
at hand) grows very quickly with the dimension iV since dimV = entier[N 2 / 4] in 
general. This makes the present N = A model quite exceptional because in the very 
next N = Q model one already has dimV = 9, etc. 

In the purely formal setting, a sufficiently well-motivated reduction of the number 
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of the "relevant" matrix elements should be proposed in the future, therefore. The 
very first steps in this direction have only been made very recently - in ref. [23] certain 
additional symmetries have been introduced via certain non-Hermitian parity-type 
operators V ^ V\ etc. 

In the more realistic considerations the relevance of the present model relates to 
the situations where some of the energy levels of a quantum system get close to each 
other. A number of experimental as well as theoretical challenges is encountered. On 
one side, during a variation of parameters the so called avoided level crossings may 
be observed in some nonrelativistic systems like atomic nuclei [24]. On the other 
side, a confluence of the two energy levels (at a point of &D) may be followed by 
their subsequent complexification. 

In the vicinity of a point of dT> a nontrivial innovation of the physics of the model 
is often needed in its phenomenological applications. For illustration we may recol- 
lect an electron in a critically strong field where the single-particle Dirac equation 
must necessarily be replaced by its field-theoretical extension including many new 
degrees of freedom [25]. In a related brief comment [26] we emphasized that even 
on the level of the practical analyses of quantum systems using some oversimplified 
phenomenological models it is not always easy to draw the clear separation line be- 
tween the avoided and unavoided level crossings. A reliable separation of the two 
seem strongly model-dependent at present. All the future extension of the scope of 
the quantitative analysis of the models will be welcome, therefore. 

A deeper study of the phenomenon of the complexification of the energies to larger 
dimensions will be well motivated not only by its purely mathematical appeal but 
also by the very pragmatic needs of a clarification of the possible and eligible patterns 
of the spectra in phenomenological models. In this sense, our present selection of the 
specific illustrative VT— symmetric Hamiltonians H in a certain "first nontrivial" 
matrix form may be perceived as a natural starting point of such an effort. 
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Figure captions 



Figure 1. The centered circle (19) and the two hyperbolas 

C(a,b,d,d) = with the respective centers at (a, 6) = (—1,3) 
and (a, b) = (1, —3) (marked by medium circles) in a — b plane 

at c 2 = d 2 = 1.6 

Figure 2. Spiked shape of the physical domain V(a,c) near its 
lower left corner 
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